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Abstract. The extremely luminous supernova SN2006gy is explained in the same way as other 
SNIIn events: light is produced by a radiative shock propagating in a dense circumstellar envelope 
formed by a previous weak explosion. The problems in the theory and observations of multiple- 
explosion SNe Iln are briefly reviewed. 
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INTRODUCTION 



The discovery of SN2006gy yj, |2|] demonstrates that some supernova (SN) events pro- 
duce 10 or even 100 times more visible photons than other, already powerful explosions. 
The anomalously high power of the emission of SN2006gy demands an explanation. 

SN2006gy is of SNIIn type and it revived interest in SN models where light is 
produced by a long living radiative shock which propagates in a dense circumstellar 
envelope. 

I discuss problems in the theory of SNIIn events and in the physics of supercritical 
radiative shocks. The powerful visible light of SN2006gy can be easily explained by a 
radiative shock bom due to a collision of SN ejecta with a dense cloud formed by a weak 
explosion some years before the SN. Strong X-ray emission of SNIIn near maximum 
light may be absent since it is absorbed by the dense cloud or not produced at all in the 
radiation-dominated shock. 



SUPERNOVA TYPES 



The models with a long living radiative shock running in a dense circumstellar envelope 
were invoked earlier P, 0] to explain the unusual properties of other powerful super- 
novae with narrow emission lines in their spectra, SNIIn. SN2006gy also belongs to the 
same SNIIn class. 

A simple diagram below illustrates the relation between different astronomical types 
of supernovae which are classified purely on the appearance of their spectra near maxi- 
mum light, irrespective of the underlying physics. For example, we believe that SNe II 
are born when a giant star with an H-rich atmosphere has a powerful explosion in its 
core. This explosion may be a result of a catastrophic collapse of the stellar core. The 
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details of the explosion of core-collapsing SNe are still not clear, in spite of many im- 
portant results obtained by workers in this field. 




Hydrogen No Hydrogen 

Narrow Lines, type Iln type la, Ib/c 

Taking into account the physics of SN explosion we can produce a couple of other, a bit 
different, diagrams for SN types. I want to subdivide the SN physics into two groups of 
problems. First, here is a diagram for the mechanism of explosion. 



Core Collapse r.type II, Ib/c 



Thermonuclear Thermonuclear 

Nondegenerate, Pair Instability, type Iln Degenerate, type la 

And the next diagram illustrates main ways to produce light during supernova explo- 
sions. Here S is entropy, and 'heating' is entropy production. 



Cooling 
no source of S 




Permanent heating Fading heating 

by shock, type Iln by radioactivity, type la, Ib/c 

Remarkably, the details of the explosion are not important for explaining the light curves 
of many supernovae that are bom from giant stars that retain their huge hydrogen 
envelopes. Successful SNII light models were already constructed four decades ago 
(largely by the work of the Soviet group in Moscow and Riga [5]) due to this insensitivity 
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TABLE 1. Four kinds of deaths for non-rotating stars 



Main Seq. Mass 


He Core 


Supernova Mechanism 


10 < M < 95 


2 < M < 40 


Fe core collapse to a neutron star or a black hole 


95 <M< 130 


40 < M < 60 


Pulsational pair instability followed by Fe core collapse 


130<M<260 


60 <M < 137 


Pair instability supernova 


260 <M 


137 <M 


Black hole. Possible GRB ? 



to details. The light of SNe II is the manifestation of entropy produced during a short 
period (hours to days) of the shock propagation in the body of the presupernova star. 
We cannot exclude the possibility that in some rare cases the SNe II are produced by 
thermonuclear explosions, not by a collapse, inside hydrogen envelopes. This may be 
the case for SNe Iln, and especially, SN2006gy. A supernova of type II shines for several 
months thanks to the heat stored in its body (the shock dies quickly), while in an SNIIn 
the heat (entropy) is replenished by the shock living several months. 

If the hydrogen is lost by a massive star, then we have an SNIb/c. The same (still 
unknown in details) core-collapse mechanism may lead to their explosions, as in SNe 11, 
however, the light is due now to entropy produced by radioactivity: the decays ^^Ni 
Co — Fe which lead to a slower heating of ejecta. This way of producing light is 
most important for SNe I of all subtypes. However, some contribution of radioactivity 
is clearly present in late light curves of type II supernovae as well, and for SN1987A in 
LMC it was important already before its maximum light, a month after the explosion. 

If the radioactive mechanism was responsible for the light of SN2006gy, then the 
amount of ^^Ni must be higher than IOMq [6, 7, 8] . This immediately implies a very 
large mass for the presupernova star, more than lOOM©. More important, this implies 
a huge explosion energy, (50 — 80) x 10^^ erg. One can delineate four kinds of deaths 
for massive stars, see Table \T\ There is some uncertainty about the exact values due 
to uncertainties in mass-loss, rotation etc. Such massive stars do exist, and they can 
have powerful explosions in their oxygen cores which experience instability due to the 
creation of large numbers of electron-positron pairs. 



PAIR INSTABILITY SUPERNOVAE 

The word 'instability' refers here to hydrodynamics, to mechanical equilibrium, not to 
the process of pair creation which is quite stable and reversible in a thermodynamic 
sense in stellar interiors. A massive star loses its mechanical stability when pairs are 
being created because the adiabatic exponent 7 goes down at T ~ 0.1 MeV , see Fig.[I} 
the work of contraction is spent in creating new particles and not for raising the momenta 
of particles that already exist and which provide for the equilibrium pressure. 

We can easily estimate the path that leads the star into the domain of the pair-creation 
instability. From the virial theorem, for a star of mass M and radius R, omitting all 
coefficients of order unity, PcV ~ PcR^ ~ Gj^M^/R. Hence, the pressure Pc in the center is 

Pc ~ G^M^/R^, while the density pc'^M/R^, and they are related as Pc ~ G-^M'^I^ pc^^ . 
So, if we have a classical ideal plasma with P = MpT / pi, where M is the universal gas 
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FIGURE 1. Left: adiabatic exponent in the low density asymptotic limit taking into account pair 
creation Jgl], see also ifToll . Right: evolution track for one of the models ifTlll with initial M ~ lOSM© 
(solid line). The approximate boundary of pair-instability is shown by the dashed line 



constant, and /i - mean molecular mass, we find 

If we have a relativistic addition of aT^ to P, the law Tc °^ pl^^ is the same (but the 
coefficient is a bit different). These relations have been already shown at this conference 
by Naoki Yoshida and by Marco Limongi. The right panel of Fig.[T]shows how a massive 

1/3 

Star follows the law Tc °^ pc until entering the domain of pair- instability. We had 
already computed the light curves of pair- instability supernovae for M > 130 Mq (third 
line in Table [1]) some years ago (with S.Woosley and A.Heger), but we do not like them 
in the case of SN2006gy because we do not see the evidence for a tremendously high 
explosion energy in that case. If the explosion energy were 2 orders of magnitude higher 
than for normal SNe, then it should be seen in very broad spectral lines. What we see in 
SN2006gy is different: it has narrow P Cyg lines (hence, it is type Iln) superimposed on 
moderately broad emission component (~ 5 x 10^ km/s [2]). There is no sign of huge 
kinetic energy in this event. 

Multiple ejections in SNIIn 

Supernovae of type Iln are among the most powerful transients in visible light. No 
radioactive material is needed to explain their light during the first several months: the 
light is produced by a long living radiative shock in a dense circumstellar envelope. This 
is the main difference with standard SNe II where the shock breaks out into rarefied 
interstellar medium and disappears quickly. Spectra and light curves of SNIIn can be 
explained only when the number density of circumstellar matter at radii of ~ 10^^^^^ cm 
(where the narrow lines are formed) is unusually high, like 10^^^*^ cm^^, see Fig. [21 This 
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implies that a large mass (on the order of Mq and larger) must be ejected within years, 
or even months, before the observed SN. The first ejection may have kinetic energy 
appreciably lower than a standard supernova. The slow motion of its matter explains the 
narrow lines of a type Iln SN. The supernova itself is an explosion of a normal energy, 
but inside a dense cloud. 

Paper [3] was the first to suggest that an SNIIn had a precursor, a relatively weak 
explosion ejecting a large slowly moving mass. A dramatically high mass loss needed 
for the formation of a dense envelope shortly before SNIIn 1994W (Fig. [2l) has been 
derived in [4]. SNIIn 1995G is explained in [12] in a model similar to the one presented 
in [4]. Double explosions may be observed also for other SN types [13]. 

According to an ~ 11 Mq star might produce strong flashes in the semi- 
degenerate 0-Ne-Mg core a few years prior to SN explosion and the strongest flash 
could eject most of the hydrogen envelope with velocities ~ 100 km/s. More recent 
evolutionary computations do not support the conclusion on strong Ne flashes, but this 
complicated problem deserves further investigation. 



Pulsational pair-instability 

Another workable mechanism for multiple-explosion SNe has been proposed in our 
paper ill]] to explain SN2006gy. It works for a high initial mass of the presupernova star 
~ 1 IOM0 (see Fig.O. This mechanism is also based on the pair-creation instability, but 
there is no catastrophic collapse or full explosion of the star. 

The second line of TableiHshows that between 95 and 130 Mp) a relatively unexplored 



phenomenon of pulsational pair instability supernova il5Lll6l ll7[] occurs. An instability 
in the mechanical equilibrium is encountered, as in the heavier stars, during the evolution 
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along Tc p^y path (Fig. [!])• A thermonuclear explosion of oxygen occurs, but the 
energy released is inadequate to unbind the entire star. It suffices, however, to eject many 
solar masses of surface material, including the hydrogen envelope, in a series of giant 
'pulses' (explosions of various strength). 

The binding energy for the hydrogen envelope of (95 — 130)Mq stars is ~ (0.1 — 
1) X 10^^ erg while the energy of an explosive pulse is about two orders of magnitude 
higher. The velocity is in the range 100 5000 km/s depending on amount of explosive 
burning and mass ejected. After a pulse, the remaining core contracts searching for a 
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new equilibrium state. It obeys the °^ pc law again, but now the mass is lower, so 
the track is a bit different (see Fig.[T]). The time required for the contraction is sensitive 
to the strength of the first pulse. If it cools down severely after the pulse expansion, it 
may be centuries before the star ignites burning again. If it remains hotter than 1 .5 x 10^ 
K, it may only take days. After one, two or several explosion pulses the remnant of the 
massive star continues to live, contrary to other supemovae, and eventually it should 
collapse. 

To explain SN2006gy, we consider the evolution of a star with initial mass llOM©. 
The evolution is calculated using the Kepler code [fisl. [6ll with reduced mass loss. The 
star encounters pair-instability with a total mass of 15 Mq (a helium core of 5OM0) and 
experiences the first pulse ejecting a cloud with mass ^ 25 Mq and ^kin ~ 1 .4 x 10^*^ erg. 
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FIGURE 2. Left: SN1994W sttucture fT] prototypical for SNe Iln. Right: light curve models for 
SN2006gy. Dots - observations [2J, the last observed point is from [21], courtesy M.Tanaka. The solid 
line is the model discussed in the text (with numerical resolution higher than in fill] '), and the dashed one 
where the velocity of all the ejecta has been multiplied by 2 (hence an artificial increase in the explosion 
energy to 2.9 x 10^' erg). The dotted line is for the doubled density in ejecta 



About 7 years later the remaining core produced the second explosion with ^ IOMq and 
Ekin ~ 7.2 X lO^'^ erg. Thus, the second ejection was faster and it had sent a shock wave 
into the massive cloud. The light from this radiating shock we see as SN2006gy (Fig. [J]). 
The shock is still inside the cloud, it has not yet broken out, more than one year after the 
explosion of SN2006gy. 

SOFT X-RAY? 

There were arguments against the model of a shock moving into shells from previous 
explosions for SN2006gy. One can estimate the temperature of the shock using standard 
formulas for a rarefied medium and find that it must be very high, so SN2006gy must be 
a powerful source of X-ray emission. Chandra X-ray observatory has measured the X- 
ray flux on 2006 Nov 14 from the direction to SN2006gy. A best fit for X-ray luminosity 
in (0.5-2 keV) is 1.65 x 10^^ erg/s ^ , but it is a few orders of magnitude lower than 
expected in the naive estimates. 

There is no problem with this in our model: the cloud of 25 solar masses is almost 
transparent to the visible light of the shock because it is almost neutral, but exactly due 
to this reason it is fully opaque to X-ray light. A large mass lies above the shock, see 
Fig. [3l The zero of the mass coordinate in Fig. [3] is the inner edge of the ejecta. The 
radiative shock is located at Mr ~ 6Mq, or logr ^ 15.5. 

If we look into other data on SNe Iln we see that all of them are discovered late in 
X-rays! See Table [2l This is no wonder in our models: the shock is buried in the material 
of the first ejection. 
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FIGURE 3. Hydrodynamic structure at 70 days as a function of radius (left panel) and of a Lagrangean 
mass coordinate (right panel). Black solid line is density (left Y-axis). Violet solid line is logarithm of 
absolute value of luminosity IL40I (units 10^'' erg/s). Dotted blue line is velocity vg (units 10^ cm/s), green 
line is logarithm of temperature T (in K), and the red line is Rosseland optical depth t. The scale for 
logL4o, vg, log T, and t is on the right Y-axis. 



TABLE 2. Data on X-ray observations of 
SNe Iln Q 



SN 


galaxy 


distance 


discovery 






Mpc 


day 


1986J 


NGC 891 


9.6 


3,300 


1988Z 


H-03-28-022 


89 


2,370 


1994W 


NGC 4041 


25 


1,180 


1995N 


-2-38-017 


24 


440 


1998S 


NGC 3877 


17 


678 


2005ip 


NGC 2096 


30 


490 


2005kd 


PGC14370 


64 


450 



An alternate explanation for the low X-ray flux is possible. After the work of physi- 
cists on nuclear explosions in atmosphere in the 1940s and 1950s, we know that the 
preheating effect becomes so large in strong supercritical shocks that the viscous jump 
in pressure and density diminishes and completely disappears. In radiation dominated 
shocks not only the preheating effect is important; in addition the momentum transfer 
from photons to electrons (and hence to ions, via the electric field) is very large. This also 
destroys the viscous jump at the shock front. Imshennik and Morozov [19] have found 
with an accurate accounting of the photon transfer that this happens when Pr/Pg — 8.5. 
This effect implies that postshock temperature may be so low that it does not attain keV 
range while inside the envelope. All the heat is taken away by 'cold' photons. Thus, the 
Chandra results [2] tell us something about the interaction of the the first pulse ejecta 
with the ISM, not about the main shock penetrating the pre-ejected shell, which shines 
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in visible light as for SN2006gy. 

CONCLUSIONS 

1. Radiating shocks are the most probable sources of light in most luminous super- 
novae like SN2006gy. 

2. The medium in which the shining shock propagates is naturally produced in mas- 
sive star evolution due to violent non-linear pulsations when e e+ pairs become 
appreciable in the pressure in their interiors. 

3. The supercritical radiative shock must be well below the X-ray temperature. Even 
if the shock becomes hot, the overlaying matter absorbs X-rays for a long time. 

4. If SN2006gy is a pulsational-pair-instability SN, a better understanding of mass- 
loss is needed. Otherwise, one has to find other evolution paths to SNe with double 
(or multiple) explosions. 
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